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ABSTRACT 

For nearby K dwarfs, the broadening of the observed Main Sequence at low metal- 
licities is much narrower than expected from isochrones with the standard hehum- 
to-metal enrichment ratio AY/AZ 2. Though the latter value fits well the Main 
Sequence around solar metallicity, and agrees with independent measurements from 
HII regions as well as with theoretical stellar yields and chemical evolution models, 
a much higher AY/AZ 10 is necessary to reproduce the broadening observed for 
nearby subdwarfs. This result resembles, on a milder scale, the very high AY/AZ 
estimated from the multiple Main Sequences in wCen and NGC 2808. Although not 
"inverted" as in wCen, where the metal-rich Main Sequence is bluer than the metal- 
poor one, the broadening observed for nearby subdwarfs is much narrower than stellar 
models predict for a standard helium content. We use this empirical evidence to argue 
that a revision of lower Main Sequence stellar models, suggested from nearby stars, 
could significantly reduce the helium content inferred for the subpopulations of those 
globular clusters. A simple formula based on empirically calibrated homology relations 
is constructed, for an alternative estimate of AY/AZ in multiple main sequences. We 
find that, under the most favourable assumptions, the estimated helium content for 
the enriched populations could decrease from Y ~ 0.4 to as low as y ~ 0.3. 

Key words: stars: abundances - stars: fundamental parameters - stars: subdwarfs - 
Hertzsprung-Russell (HR) diagram - globular clusters: individual: uj Cen, NGC 2808 



1 INTRODUCTION 

Helium is the second most abundant element in the Uni- 
verse, having been produced by Big Bang Nucleosynthe- 
sis (BBN) with a universal mass fraction of Yp ~ 0.24 
(Steigman 2007); after that, successive stellar generations 
synthesize metals (Z) and an additional amount of helium 
Ay. The characteristic helium-to-metal enrichment rate 
AY/AZ is estimated to be around 2, both observationally 
from HII regions (e.g. Fukugita & Kawasaki 2006; Peimbert 
et al. 2007; Izotov, Thuan & Stasinska 2007) and theoret- 
ically, from models of stellar nucleosynthesis and chemical 
evolution of the Galactic disc (e.g. Chiosi & Matteucci 1982; 
Maeder 1992; Carigi & Peimbert 2008; Casagrande 2008) 

In spite of its large abundance, helium is an elusive ele- 
ment. It can be measured directly, from spectroscopic lines, 
only in stars hotter than ~ 10,000 K: young, massive stars 
(or their surrounding HII regions) or blue Horizontal Branch 
stars (where though, apart from a small temperature range, 

* E-mail: Iporti, cflynn@utu.fi, lucaOmpa-garching . mpg . de 



the surface abundance may not trace the original one, due to 
helium sedimentation and metal levitation: Michaud, Vau- 
clair & Vauclair 1983; Michaud, Richer & Richard 2008; 
Villanova, Piotto & Gratton 2009 and references therein). 
Only indirect methods can be used for stars of lower mass 
and cooler temperatures, which constitute the bulk of the 
Galaxy's stellar population. One such method relies on the 
broadening of the lower Main Sequence (MS): in fact the 
location of low-mass, long-lived stars in the HR diagram 
depends on their metallicity Z and their helium content Y 
— or equivalently, on AY/AZ. Concomitant increases in Z 
and Y shift the MS in opposite directions, so for a given 
variation AZ the two MSs are more spread apart (the metal 
rich being cooler at a given luminosity) if the correspond- 
ing Ay is lower; conversely, if Ay is large "enough", the 
MSs may even invert, with the more metal (and helium) 
rich being bluer (i.e. hotter; e.g. Fernandes et al. 1996). 

This method dates back to Faulkner (1967); early stud- 
ies based on ground-based parallaxes deduced a large and 
quite uncertain value of AY/AZ ~ 5 ± 3 from the apparent 
overlap of MSs of all metallicities (e.g. Perrin et al. 1977; 
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Fernandes et al. 1996). The situation greatly improved with 
Hipparcos parallaxes, so that a net separation of the MS as 
a function of Z could be detected in the HR diagram. Stud- 
ies based on Hipparcos distances concluded Ay/AZ = 2 — 3 
(Pagel & Portinari 1998; Jimenez et al. 2003). More recently, 
Casagrande et al. (2006, 2007) further improved the analysis 
by compiling a much larger sample of K dwarfs with homo- 
geneous multi-band photometry, for which bolometric mag- 
nitudes and effective temperatures were derived with the 
InfraRed Flux Method (IRFM). The interpretation of the 
broadening of the lower MS could thus be performed in the 
theoretical HR diagram (Mboi vs. log T^a) where comparison 
to stellar models is more straightforward and the effects of 
/SY/ /S.Z more prominent (Castellani, Degl'Innocenti & Mar- 
coni 1999). Around solar metallicities this analysis yielded 
Ay/AZ ~ 2; but at lower Z Casagrande et al. (2007) found 
that the observed Main Sequence is narrower than expected 
from stellar models computed under the standard assump- 
tion Ay/AZ 2, thus implying a very steep Ay/AZ ~ 10 
(Fig. [5]). This change in slope is at odds with Galactic chem- 
ical evolution models, that predict a ratio substantially con- 
stant with Z (Carigi & Peimbert 2008; Casagrande 2008); a 
constant AY/ AZ ~ 2 is also found for both metal-rich (e.g. 
Peimbert 2003; Balser 2006) and metal-poor (e.g. Fukugita 
& Kawasaki 2006; Peimbert et al. 2007; Izotoz et al. 2007) 
HH regions. More importantly, AY/ AZ ^ 10 at low metal- 
licities implies a helium content for nearby subdwarfs much 
lower than Yp, in awkward constrast with the cosmological 
floor set by standard BBN. 

A combination of small sample size and Z-dependent 
bolometric corrections in the observational HR diagram had 
masked this effect in previous studies of AY/ AZ, although 
very low Y values had already been noticed in a handful of 
low metallicity stars with available IRFM Mboi and Tcff (Fer- 
nandes et al. 1998; Lebreton et al. 1999). This result points 
toward inadequacies in Main Sequence stellar models of low 
metallicity, and we now plan to investigate the issue fur- 
ther, as this may lead to reconsideration of the problem of 
the helium enrichment in globular clusters. 

In fact, at absolute magnitudes comparable to those of 
the local stars studied in Casagrande et al. (2007, Mv ^ 
5.5 — 7.5), multiple Main Sequences have been discovered in 
some globular clusters and interpreted as evidence for huge 
helium enhancement in a sub-population (Bedin et al. 2004; 
Norris 2004; Piotto et al. 2005, 2007); see Fig. [T] The re- 
quired AY / AZ > 100 is extremely difficult to explain 
with stellar nucleosynthesis and chemical evolution models 
(Karakas et al. 2006; Prantzos & Charbonnel 2006; Maeder 
& Meynet 2006; Choi & Yi 2007, 2008; Romano et al. 2007, 
2009; Renzini 2008; Yi 2009; Marcolini et al. 2009; Peng & 
Nagai 2009). 

The helium content of the sub-population also has a ma- 
jor role in shaping other regions of the HR diagram (Catelan, 
Valcarce & Sweigart 2009a and references therein); in par- 
ticular it is reflected in the morphology of the Horizontal 
Branch (D'Antona & Caloi 2004; Caloi & D'Antona 2005, 
2007; Lee et al. 2005). While it is likely that a significant 
helium enhancement is present in those stellar populations, 
the purpose of this work is to show that the revision of 
low metallicity MS stellar models, needed to cure the prob- 
lem of the high AY / AZ and sub-primordial Y deduced in 
nearby K dwarfs, may significantly reduce current estimates 
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Figure 1. Le/t panel: HR diagram of nearby K dwarfs (filled 
circles) from the sample of Casagrande et al. (2007), overplotted 
on the MSs photometry of ojCen (SoUima et al. 2007), adopting 
E{B-V) = 0.11 (Lub 2002) and dereddened distance modulus of 
13.7 (Bellazzini et al. 2004; Del Principe et al. 2006). The maxi- 
mum split between the red and blue MS of ujCen A{B — R) = 0.1 
occurs around Mn ~ MboI ~ 6.5 (SoUima et al. 2007). Right 
panel: nearby K dwarfs in the theoretical plane with isochrones 
having Z = 0.001,^0,0.04 and AY/AZ = 2 overplotted. Clear 
mismatch in the broadening appears for low metallicities. Below 
A^Bol = 5.4 the broadening of the MS is age independent. We 
also show the reddening/extinction correction corresponding to 
E{B — V) = 0.10 that would be needed to reconcile the metal- 
poor stars with theoretical isochrones. 



of AY/AZ in globular clusters with multiple MSs, easing 
their theoretical interpretation. 

We show this by means of an exercise based on ho- 
mology relations. Theoretical homology relations, describ- 
ing how the location of the low MS in the HR diagram 
depends on metallicity and helium content, are briefly re- 
called in Section [21 In Section 3, we apply homology relations 
to interpret the split of multiple MSs in globular clusters, 
and show that the inferred helium enrichment for the blue 
population(s) agrees with the results of detailed isochrone 
analysis. In Section 4 we confirm that, in general, homology 
relations render very well the behaviour of stellar models as 
a function of the helium content. However, both isochrones 
and homology relations fail the interpretation of the low 
MS of nearby stars — as we cannot accept sub-primordial 
helium contents for low metallicity stars. Therefore, in Sec- 
tion 5 we proceed to calibrate empirical homology relations, 
that return the expected AY/AZ = 2 for the low MS of 
nearby stars with Z < Zq. Such empirically calibrated ho- 
mology relations also return, for the helium-rich subpopula- 
tions of globular clusters, a helium content Y 0.3, rather 
than the much higher Y ~ 0.35 — 0.4 of standard analy- 
sis. This lower value of helium enrichment can be recon- 
ciled far more easily with present chemical evolution mod- 
els (Karakas et al. 2006; Renzini 2008; Yi 2009). Finally, we 
conclude in Section 6 with a discussion of the possible phys- 
ical processes that could improve stellar models for low MS 
stars, so as to solve the "helium problem" both for nearby 
stars and for globular clusters. 
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Figure 2. Broadening of the ZAMS in log T^s as a function of metallicity, as predicted by homology relations (Eq. [JJ- The broadening 
is referred to the primordial ZAMS {Zp, Yp\ left panel) or to the solar ZAMS (Zq = 0.017, Yq = 0.263; right panel) for ^ ranging 
from to 10 as indicated. Dotted lines show the approximate relation proposed in Eq. |3] The thick red line highlights the case = 5, 
which corresponds to MS inversion around Zq. At increasing AY/llZ, MS inversion occurs at lower and lower Z. 



2 HOMOLOGY RELATIONS 

Homology relations (Cox & Giuli 1968), holding for radiative 
structures like the energy-producing cores of MS stars of 
M < 1 Mq , express the difference in MboI of a Zero Age 
Main Sequence of composition {Y, Z) with respect to another 
reference ZAMS of composition {Yr,Zr). In the hypothesis 
that the compositions are related by: 

AY 

Y^Yr+—{Z-Zr) 

the formula is: 
AMboi= -1.59 log [l- 

-3-33 log [l- (3+5x:iz.) (^-^0] 

r 1 
-0.8675 log ^l^^^{Z^Zr)\ 

-0.8675 log (^12-±i_) 

where (5=1 + and Xr — 1 — Y,. ~ Zr- Notice that the 
same relation holds also for a generic combination of {Y, Z) 
replacing S (Z - Zr) {Y - Yr) + {Z- Zr). The first term is 
related to the rate of thermonuclear energy generation, the 
second term to the molecular weight, and the final two terms 
to opacity (Cox & Giuli 1968; Fernandes et al. 1996; see also 
the Appendix). For Zr = Q, Yr = Yp (primordial helium 
floor due to BBN) this formula corresponds to eq. 3 of Pagel 
& Portinari (1998). Another typical approach is to refer to 
the solar ZAMS, as the solar model is the basic calibration 
of stellar tracks and isochrones; then {Yr, Zr) = {YQ, Zq). 

In the case of old stellar populations, because evolution 
affects luminosities approximately brighter than Mboi=5.4, 



it is more useful to translate the broadening in Mboi into 
a broadening in Tcff using the slope of the lower MS. In 
the Padova isochrones used in Casagrande et al. (2007), this 
slope is about 17 mag per dex in log(Tcff), henceforth: 

A log Teff = -0.0935 log [l - ^ (^ - Z,)] 

-0.196 log fl - ^^^l±t^(Z - Zr)] 

' (2) 

-0.051 log ^1-^^{Z-Zr)\ 

-0.051 log 

(cf. eq. 4 of Pagel & Portinari, with slightly different coeffi- 
cients due to the different adopted slope of the MS). Fig. [5] 
shows the broadening of the MS as a function of metallic- 
ity, as predicted by homology relations. For low AY/AZ 
(e.g. AY/AZ = 0) the MS is cooler at increasing metal- 
licity, as "normally" expected. When AY/AZ = 5, MS in- 
version is expected around solar metallicity (i.e. overlapping 
MS, as was deduced from pre-Hipparcos data). At increas- 
ing AY/AZ, MS inversion — namely, more metal-rich MS 
being hotter due to the overwhelming effect of the helium 
excess AY — occurs at lower and lower Z. Indeed, extremely 
high AY/AZ are necessary to interpret the MS inversion at 
the low Z oi uj Centauri. 



3 HOMOLOGY RELATIONS AND GLOBULAR 
CLUSTERS 

Evidence for multiple stellar populations in some Globular 
Clusters comes from turn-off, subgiant and red giant branch 



4 Portinari, Casagrande and Flynn 



splits; bi- or multi-modal distribution of light elements; and 
the presence of multiple, distinct MSs (e.g. Piotto 2009). In 
the latter case, because of the faintness of the stars, deep 
photometric observations are needed and clear evidence for 
multiple MSs is currently limited to tjCen and NGC 2808; 
in addition, evidence for MS broadening has been recently 
found in 47 Tuc (Anderson et al. 2009). 

In this section we analyze the multiple MSs of the 
above mentioned clusters by means of homology relations, 
and compare the results to those obtained with modern 
isochrone fitting. 

3.1 u Cen 

Detailed isochrone analysis indicates that a helium enrich- 
ment of Ay ~ 0.15 between the two populations is needed 
to reproduce the inverted MSs of tjCen (Norris 2004; Pi- 
otto et al. 2005; Lee et al. 2005; SoUima et al. 2007). We 
show here that homology relations yield a very similar re- 
sult. To apply homology relations (Eq. [2]|, we must trans- 
late into a split in log T^s the maximum observed colour 
split between the blue (bMS) and red (rMS) main sequence: 
A(i3 ~ R) — 0.1, which occurs at a dereddened B — R — 1.2 
for the rMS (SoUima et al. 2007). We apply the colour- 
temperature-metallicity scale of Casagrande et al. (2006) 
that was derived for nearby K dwarfs, adopting [Fe/H]j,^g = 
-1.6 and [Fe/H]j^[^g = -1.3 (SoUima et al. 2007). We ob- 
tain for the rMS rcff=5120 K and A log Tea = 0.0185 be- 
tween the two sequences. Considering possible uncertainties 
in the reddening which affect the absolute B — R values 
(while the differential A(_B — R) — 0.1 is robust) we esti- 
mate Alog(reff) = 0.0185 ± 0.0015. Adopting the updated 
and extended temperature scale by Casagrande et al. (2010) 
also yields AlogTcfi estimates within this range. 

To derive the metal mass fraction Z that is relevant 
for homology relations, we adopt for both sequences an en- 
richment [a/Fe] = -1-0.3 as supported by various spectro- 
scopic studies (SoUima et al. 2007 and references therein; 
Villanova et al. 2007). Using the relation of Yi et al. (2001) 
to compute the global metallicity leads to [M/H]^j^g = 
-1.37 and [M/H]^^jg = -1.07. Simple scaling with Zq 
returns Z^ms = 0.00072 and Z^ms = 0.00144. We fur- 
ther assume Y^ms = 0.246 — but note that the resulting 
Ay = ybMS — is not sensitive to the adopted value 

of yMS. Fig. [3] shows the split in AlogTcfi predicted by 
homology relations with respect to the rMS, i.e. applying 
Eq. [2] with (yr, Zr) = (yMS, ^rMs)- The observed split be- 
tween the rMS and bMS is reproduced for AY/AZ ~ 200, 
implying Ay = 0.144 and ybMS = 0.39. These values for 
Ay and ybMS are in excellent agreement with the results ob- 
tained by Piotto et al. (2005), Lee et al. (2005) and SoUima 
et al. (2007) by means of detailed isochrone analysis — while 
the value for AY/AZ may vary significantly, depending on 
the assumed difference in metallicity AZ in various studies. 
For instance, by adopting Z^ms = 0.001 and ZbMS ~ 0.002 
as in Piotto et al. (2005), homology relations reproduce the 
observed temperature split with AY/AZ ~ 150, implying 
again Ay ~ 0.15 as a robust result. A rigorous transforma- 
tion from [M/H] to Z (e.g. Casagrande et al. 2007) should 
take into account also the very different helium fraction 
in the two sequences (yMS = 0.246 and ybMS = 0.39), 
returning Zrus = 0.00076 and ZbMS = 0.00120; homol- 
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Figure 3. Temperature split of the MS predicted by homology 
relations compared to the split between the rMS (dot, taken as ref- 
erence MS in Eq.[21l and the bMS (star symbol) of ojCen. The ob- 
served temperature split is Alog(T(,fi) = 0.0185 ± 0.0015, and the 
assumed metallicities are ^rMS = 0.00072 and .ZbMS = 0.00144 
(see text). Lines indicate homology predictions for Ay/AZ in- 
creasing from to 200 in steps of 20. The observed split between 
the rMS and bMS is reproduced for AY/AZ ~ 200, correspond- 
ing to Ay = 0.144 and YbMS = 0.39. 

ogy relations then yield AY/AZ ~ 320, corresponding to 
Ay = 0.14. Thus, Ay = YbMS - y-MS is a robust resuU of 
the analysis. 

It is worth underlining that, in any case, in tjCen 
AY/AZ >> 70, this value being correctly indicated by Pi- 
otto et al. (2005) as a lower limit, but often quoted in the 
literature as the preferred AY/AZ. Values of AY/AZ 200 
or higher are implied by the current metallicity measure- 
ments of the rMS and bMS (as also in SoUima et al. 2007), 
which renders the theoretical interpretation even more trou- 
blesome (see e.g. Yi 2009). 

3.2 NGC 2808 

Another striking example of helium enriched multiple Main 
Sequences is NGC 2808, where three MSs are found at vir- 
tually the same metallicity [Fe/H] = —1.1 ± 0.03 dex (Car- 
retta et al. 2006; see also Salaris et al. 2006) and helium 
abundances of y = 0.248 (the value assumed for the red 
MS, which includes the bulk of the population), Y — 0.30 
(middle MS) and Y = 0.37 (blue MS; D'Antona et al. 2005; 
Lee et al. 2005; Piotto et al. 2007). In this case, formally 
AY/AZ ^> oo (or, AY/AZ > 800) and one should just 
focus on Ay. Also in this case, we show that homology re- 
lations yield an estimate of AY in good agreement with 
isochrone analysis. 

As for wCen, the scale of Casagrande et al. (2006) was 
used to convert the colour split into a temperature split. In 
NGC2808, the maximum split of the bMS and rMS with re- 



/^Y / /S.Z in Globular Clusters: insight from nearby stars 5 



0.02 



0.01 - 





1 1 1 1 1 1 1 1 




^ — 


— — 




"~ ■ ~ — ■ AXsOJ. 

^^H~~ 






AY^^O ^ 








0,001 



0,0015 



0.002 
Z 



0.0025 



0.003 



Figure 4. Temperature split of tiie MS predicted by liomology 
relations compared to the split between the red MS (dot, taken 
as reference MS in Eq. [2J, the middle MS (square symbol) and 
the blue MS (star symbol) of NGC 2808. In this figure homology 
relations are marked as a function of Ay rather than AY/AZ 
(see text), with Ay increasing from to 0.17 in steps of 0.01. 
The observed splits are reproduced for AY (mMS-rMS)=0.074 
and Ay(bMS-rMS)=0.143, corresponding to Y^ms = 0.322 and 
^bMS = 0.391. Errors bars are estimated converting in AlogT^ff 
an uncertainty of 0.02 magnitudes for the value of A{B — I). 



spect to the mMS (±0.08 mag) occurs around a magnitude 
?TiF8i4w = 21.3, where the colour of the mMS is F475W — 
F814W = 1.72. Correcting for reddening (E(B-V)=0.18, Pi- 
otto et al. 2007) and transforming the WFC/ACS Vegamag 
system photometry into Johnson-Cousins according to Siri- 
anni et al. (2005), these colours and broadening correspond 
to B - / = 1.74 and A{B - I) ± 0.11 mag, equivalent to 
Te// = 5020 K for the mid MS and AlogTeff = ±0.01. 

As AZ vanishes and AY/ AZ diverges in the case of 
NGC 2808, it is convenient to recast the homology relation 
in Eq. [2] replacing 5 {Z - Zr) ^ {Y - Yr) + {Z - Zr). Fig. g] 
compares the predicted temperature split as a function of 
Ay to the observed one. For an assumed Y-^ms = 0.248, 
the resulting FrnMS = 0.32 and YbMS = 0.39 are again in 
excellent agreement with isochrone analysis (see Fig. 2 of 
Piotto et al. 2007). 

3.3 47 Tuc 

Theoretical homology relations agree very closely with de- 
tailed isochrone analysis in the estimate of the Ay and 
AY/ AZ, characterizing the split of the multiple MSs of 
uoCen and NGC 2808. We have checked that this holds true 
also for the (less extreme) case of 47 Tuc. The observed 
broadening of the MS of 47 Tuc is A(F606W" - F814W") = 
0.013 around an average colour F606VK - F814W^ ~ 0.9 
at F&Q&W = 20.5; if interpreted in terms of a dispersion 
in helium abundance, this implies Ay = 0.026 (Anderson 



et al. 2009). To analyze the same dispersion with homology 
relations, we have translated the observed broadening into 
Johnson colours, and then into effective temperatures. 

With the reddening and distance modulus listed in the 
catalogue by Harris (1996; in its 2003 online version) the 
broadening amounts to AiV — I) = 0.017 around Mv ~ 7.3 
and V-I = 1.04. With a metallicity [Fe/H] = -0.76 (Harris 
2003; Carretta et al. 2009), this average V — I colour corre- 
sponds to Tcft ~ 4800 K and the broadening to A log T^s — 
0.003. Homology relations (Eq. [2]) reproduce such broaden- 
ing with Ay — 0.023, again very close to the conclusions 
from isochrone fitting. 



4 HOMOLOGY RELATIONS AND 
ISOCHRONES 

We have shown that theoretical homology relations, applied 
to the multiple (or broadened) MSs of globular clusters, pro- 
vide a similar interpretation of the data in terms of Ay and 
AY/AZ, as detailed isochrone analysis. We now compare 
directly homology relations to current stellar models. We fit 
the Padova isochrones as a function of Z and AY/ AZ with 
a homology-like relation and compare the fitted relation to 
the theoretical one. 

An interpolation formula for isochrones directly in- 
spired by Eq. [2] involves at least 4 independent parameters: 
the coefficients of the first 3 terms in Eq. (2] and the coeffi- 
cient of Z within the argument of the logarithm in the 4*'' 
term. For fitting purposes, however, we found that a sim- 
plified form of the homology relations is more convenient. 
We derive it by condensing together the first three terms 
in Eq. (2] containing the dependence of the broadening on 
AY/AZ. First of all we notice that those three terms have 
a similar dependence: 



oc log 



-iz-z,.] 



{a + Xr) 

as in the second term we can approximate 
55+1 ,^ ^ , 55 



a = 0, 0.6 or 1 



(Z-Zr 



{Z - Zr) 



(3±5X, -ZO ' (3±5X,; 

Secondly, it is easy to verify that the second term of Eq. [2] 
is the leading one, contributing about as much as the sum 
of the first and third term together. Guided by these con- 
siderations, we find that the following formula: 



AlogT.ff 



-0.414 log [ 
-0.051 log 



{0.(i + Xr) 

lOOZ+1 



lOOZr 



(3) 



is a very good approximation of the rigorous homology rela- 
tions, especially for metallicities Z ^ Zq which are of inter- 
est for our discussion (Fig. [2]). We have also verified that this 
is a fully adequate approximation of the full homology rela- 
tions even for extremely high AY/AZ > 100, which are of 
interest for globular clusters, at least within narrow metal- 
licity ranges (like those separating the sub-populations of 
e.g. ujCen). This simplified formulation has the advantage 
of separating one term sensitive to the helium content, from 
the second term which depends only on metallicity. This 
clearcut separation will prove to be handy for the empirical 
re-calibration of homology relations in Section 5. 
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Figure 5. Metal versus helium mass fraction for nearby field dwarfs. Lines with two values of helium-to-metal enrichment ratio (break 
at Z = 0.015) are plotted to guide the eye. Left panel: using the isochrone fitting procedure described in Casagrande et al. (2007). Right 
panel: using numerical homology relations (see Section |4]l. 



Therefore we seek a three-parameter fitting formula for 
isochrones (and, later below, for real stars) of the kind: 

A log Teff = -A log [l - ^^^4^ [Z - Zr)] 

(4) 

-P2log(^) 

where Pi are 3 free fitting parameters. It is our experi- 
ence that this 3-parameter formula provides an adequate 
fit to isochrones, as good as (and more robust than) other 
homology-like fitting formulae with 4 or more parameters. 
We favour this form of the homology relation over the lin- 
earized one adopted by Pagel & Portinari (1998; their eq. 5) 
as it is better suited to handle large values of AY/AZ (see 
the Appendix). 

We consider both the Padova isochrones computed 
specifically for the analysis in Casagrande et al. (2007), and 
the more recent release by Bertelli et al. (2008) with vary- 
ing AY/AZ and 0.0001 Z 0.04, which fully brack- 
ets the metallicity range relevant for the present study. In 
Casagrande et al. (2007) we had checked that other sets of 
isochrones (Yonsei-Yale, Teramo, McDonald) were very sim- 
ilar in the low MS. 

We tried the fit both with respect to the solar [Zq = 
0.017) and to a very metal-poor reference isochrone {Z = 
0.0002). We explored values of AY / AZ up to 1000, ob- 
tained both with sub-solar and sub-primordial helium con- 
tents (isochrones from Casagrande et al. 2007) and with 
helium-enhanced isochrones (Bertelli et al. 2008). Very sim- 
ilar parameters are obtained in both cases, as one would 
expect from the discussion in Section [51 

We find that the behaviour of isochrones, at least in 
the range Mboi = 5.4 — 7.0 which is where the effect of 
AY/ AZ is expected to be maximal and dominant, is well 
described by homology-like relations with the following set 
of parameters: 

Pi = 0.50 ± 0.03 



P2 = 0.064 ± 0.005 
Pz = 670 ± 200 

where the values are the averages from the fits obtained 
at three fixed A/bo1 = 5.4,6.0,6.5 and Z = 0.0002 and 
0.017. In the following of the discussion, when using these 
parameters we will speak of "numerical homology relations" . 
With these parameters, the results obtained for field stars by 
Casagrande et al. (2007) are well reproduced (see Fig[5l). In 
Casagrande et al. (2007) the metal and helium mass fraction 
of the stars were computed iteratively, with a star-by-star 
isochrone fitting, whereas now Z is scaled with Zq and Y 
computed via Eq. [4] Despite the less sophisticated approach 
provided by homology relations, the overall agreement is 
good and the change in slope above and below Z ~ 0.015 is 
still reproduced. This confirms the validity of our simplified 
approach to study the broadening of the low MS. 

We notice that the fitting parameters are quite close to 
those of the (simplified) homology relation in Eq. [31 with 
the exception of P3 which is about 7 times larger. This was 
also noticed by Pagel & Portinari (1998); the corresponding 
parameter in their formalism was P3 = 1/Zo, with Zq = 0.01 
from the theoretical homology relations and Zo — 0.0015 
from isochrone fitting. The discrepancy was mostly imputed 
to a difference between the observational and the theoretical 
HR diagram (Castellani et al. 1999) while here we confirm 
it also in the purely theoretical plane. 

The parameter P3 describes the metallicity dependence 
of the opacity coefficient kq (see the Appendix) so it is no 
surprise that it may differ between detailed stellar structure 
models and simplified descriptions of kq. Opacity through- 
out the star (at least for largely radiative structures like low 
mass MS stars, where the convective envelope is extended 
in radius but not so prominent in mass, and most impor- 
tantly, energy production occurs in radiative regions) is the 
main factor determining its luminosity and henceforth its 
whole structure; so a significant difference in the parameter 
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P3 implies a significant difference in structure (or, better 
to say, in its metallicity dependence) between actual stellar 
models and homologous stars. However, P3 is unrelated to Y 
or Ay/AZ and, since the other parameters (Pi and P2) are 
similar to the expected values, in essence homology relations 
render very well the dependence of theoretical isochrones on 
variations of the helium content — as our experiment on 
wCen and NGC 2808 had already suggested. 



5 EMPIRICAL HOMOLOGY RELATIONS 

In the previous sections we have shown that theoretical 
homology relations agree very well with the isochrone- 
based interpretation of the MS splits observed in wCen 
and NGC 2808, and adequately reproduce the behaviour 
of isochrones in general, as a function of AY/A.Z. However, 
isochrones are known to fail the interpretation of the HR 
diagram of nearby metal-poor low MS stars (e.g. Lebre- 
ton et al. 1999; Torres et al. 2002; Casagrande et al. 2007; 
Boyajian et al. 2008). Therefore, in this section we define 
"empirical" homology relations, calibrated to reproduce the 
broadening of the local MS, and extrapolate the conse- 
quences for the interpretation of globular clusters. 

Casagrande et al. (2007) have shown that, while for Z > 
0.015 the broadening of the main sequence is well reproduced 
with the standard AY/AZ ~ 2, with current isochrones a 
much higher AY/ AZ ~ 10 is needed to fit lower metallicities 
(Fig. While in principle AY/AZ may not necessarily be 
constant. Fig. [5l taken at face value, also implies a helium 
content in local metal poor stars as low as Y — 0.1. Such a 
striking contrast with BBN is to be ascribed to inadequacies 
of low metallicity stellar models. 

While awaiting for a solution to this problem by 
improved stellar physics (see our final discussion), here 
we adopt a very pragmatic approach: we assume that 
AY/AZ ~ 2 also for any Z < Zq, and empirically cali- 
brate homology relations so that this result is returned for 
nearby stars. Our assumption is very reasonable, since HIT 
region measurements and chemical evolution models support 
a constant AY/AZ ~ 2 down to very low Z (e.g. Peimbert 
et al. 2007; Carigi & Peimbert 2008), as does the follow- 
ing simple argument: taking the solar bulk abundances (As- 
plund et al. 2009) and the primordial Yp — 0.240 (Steigman 
2007), one derives AY/AZ = 2.1 for Z < Zq. 

To define our empirical homology relations, let's first 
inspect the form of the relations (Eq. The first term 
includes the effects of the helium content, while the sec- 
ond term (stemming essentially from the metallicity depen- 
dence of opacity) is independent of AY/AZ. The second 
term is therefore irrelevant for the multiple MSs in GCs, 
where metallicity differences are small or vanishing. There- 
fore, if we calibrate our empirical homology relations onto 
nearby stars acting only on the first term (via the parame- 
ter Pi), we maximize the change in the role of helium, and 
hence maximize the consequences for the interpretation of 
GC'sQ 

We shall therefore start by discussing how the first term 
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Figure 6. Metal versus helium mass fraction for the dwarfs in 
the right panel of Figure [5] (zooming on metallicities around and 
below solar). Open circles are using the numerical homology re- 
lations established in Section [2] from isochrone fitting, while full 
circles are obtained using Pi = 1.5 as discussed in the text. A 
line with constant AY/AZ = 2 is also shown for comparison. 



of Eq. |4] must change, to yield an acceptable AY/AZ ~ 
2. This correspond to imputing the erroneous theoretical 
broadening of the main sequence entirely to the response 
of stellar structure to the helium-to-metal enrichment ratio 
AY/AZ. 



5.1 Empirical calibration of Pi and globular 
clusters 

Extant isochrones, or the equivalent "numerical" homology 
relations (Eq.|4]with Pi = 0.50, P2 = 0.064 and P3 = 670) 
yield, for the observed broadening A log Toff at low Z, a re- 
sulting (5=1 + AY/AZ ~ 11. As we aim at imposing 5 = 3, 
we can expect that we will need to increase Pi by a factor of 
the order of 11/3, i.e. Po = 0.50 — 2. This is easily seen 
e.g. with a Taylor expansion of the logarithm, so that the 
first term approximately goes as ~ Pi - 



(0.6-i-x^) (Z-Zr): for a 
given A\ogTcB{Z), Pi and S are inversely proportional, and 
an increase of the former by 11/3 corresponds to a reduction 
of the latter by the same amount. 

More rigorously, we keep P2 = 0.064 and P3 = 670 
fixed and optimize Pi so that the inferred helium content 
of nearby K dwarfs with Z < Zq follows AY/AZ = 2. 
This optimization yields Pi = 1.5 (similar to our simple 
expectations above) and the helium abundances shown in 
Fig.H 

We now re-interpret the split of the multiple MS of 



^ Ideally, one would use the data to calibrate all of the three 
parameters at the same time, or even to calibrate more complex 
relations such as Eq. (2] which would yield some handle on the 



corresponding physical ingredients. Unfortunately, the data is too 
noisy to allow for such a detailed calibration. 
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Figure 7. Left (right) panel: same as Figure [3] Q for tjCen (NGC 2808), but using Eq.gjwith Pi = 1.5. 



globular clusters by means of these empirical homology re- 
lations, with a re-calibrated Pi = 1.5. Fig. [7] shows that 
the blue sequence of ojCen is now fitted with Ay/AZ = 75 
and YbMS = 0.30. For NGC 2808, the middle and blue MS 
are now fitted with hY = 0.02 and 0.04 respectively, corre- 
sponding to YmMS = 0.27 and IbMS = 0.29. These new val- 
ues of the helium enrichment are significantly lower than the 
earlier estimates Y ~= 0.4 and are within reasonable reach 
of extant theories on stellar nucleosynthesis and chemical 
evolution (e.g. Yi 2009). 

To summarize: if we assume that the unique culprit of 
the erroneous broadening of the theoretical low-Z MS is the 
response of stellar structure to AV/AZ, and consequently 
recalibrate the first term in the homology relations, the es- 
timated Ay for the subpopulations of GC's is drastically 
reduced. The new estimates of the helium content of the 
subpopulations are no longer extreme; rather, they are close 
to the solar helium content or slightly larger {Y ^ 0.3) which 
makes them relatively easy to reconcile with chemical evo- 
lution models0 

5.2 Empirical calibration of the second term 

An alternative way to bring the theoretical broadening of the 
low MS in agreement with observations, is to act on the sec- 

^ This new estimate was obtained by recalibrating Pi on the 
low— Z range alone, i.e. using observed K dwarfs with Z < Zq 
(Fig. (6)1, as for Z > Zq isochrones yield the correct AY//\Z ~ 2 
as they are. If we rather want to calibrate the homology relations 
over the whole metallicity range of the sample, so that the same 
homology formula yields Ay/AZ ~ 2 between < Z < 0.04, we 
get a lower Pi = 1.1 as the optimized value. The interpretation of 
GCs remain similar: for ojCen, Ay/AZ = 100 and ybMS = 0.32; 
for NGC 2808, y^MS = 0.28 and ytMS = 0.30. 



ond term of the homology relations in Eq. 3] As anticipated 
above, this term only depends on metallicity and therefore 
its alteration will have no impact on the interpretation of 
the multiple MS of globular clusters, where metallicity dif- 
ferences are minimal or vanishing. But for nearby K dwarfs 
and subdwarfs, the metallicity range is significant and ex- 
ploring the effects of this second term is worthwhile. 

It is clear from (Eq. |4]) that parameters P2 and P3 are 
degenerate, in the sense that a given change in the second 
term can be obtained by modifying either of the two pa- 
rameters. Fig. [8] shows that, fixing Pi, the best combined 
solution for P2 and P3 lies roughly along a hyperbola. 

Therefore, we choose to discuss the role of the second 
term by optimizing P3, which is an interesting parameter as 
it is significantly different between the theoretical and the 
numerical homology relations (Section [4]). 

Keeping Pi and P2 fixed to the values of Section |4l 
we find that P3 ~ 150 is the optimal value to obtain 
Ay/AZ = 2 for nearby low-Z stars (Fig. [9]). However, while 
we do obtain /SY / /S.Z = 2 on average, compared to the case 
in Fig. [B] there is now considerable scatter and many stars 
remain with uncomfortably low helium abundances. It is in- 
triguing, though, that the optimized value of P3 is quite close 
to the theoretical value of 100. 

Notice that the cause of the discrepancy is entirely as- 
cribed to some sort of metallicity dependence now. For in- 
stance, the variation of the mixing-length with Z suggested 
in Casagrande et al. (2007) to avoid low helium abundances, 
in the framework of the homology relations can be formally 
described by the second term. The same can be said for 
any change in stellar models that would only respond to 
metallicity differences and not to differences in the helium 
content: it would have no impact on the interpretation of 
the multiple main sequences of Globular Clusters. 
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Figure 8. Left and central panel: reduced obtained fitting the K dwarfs in Casagrande et al. (2007) to Eq. |4] imposing /1Y/A.Z = 2 
and optimizing only one free parameter at a time. Right panel: contour lines for the reduced (values indicated by the labels) when 
only Pi is fixed: the best solution for (P2, P2) lies along an hyperbola. 
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Figure 9. Same as Figure|6]but when operating on P3 = 150 to 
obtain the Ay/AZ = 2. 



6 SUMMARY AND DISCUSSION 

In this paper we draw attention on the possible connection 
between the HR diagram of globular clusters with multiple 
Main Sequences (wCen and NGC 2808) and the broaden- 
ing of the low MS defined by nearby stars. Although not 
"inverted" as in oiCen, where the metal-rich Main Sequence 
is bluer than the metal-poor one, the broadening defined 
by nearby subdwarfs is much narrower than expected for a 
standard helium content and helium-to-metal enrichment 
law Ay/AZ ~ 2. At low metallicities, Ay/AZ ~ 10 is for- 



mally necessary to reproduce the observed broadening, but 
this is unacceptable as it implies helium fractions much lower 
than the cosmological BBN floor. 

ft is worth remarking here that the results by 
Casagrande et al. (2007) on the broadening of the low MS 
were corroborated by the analysis of a number of binary 
stars, whose masses were well reproduced by theoretical 
isochrones — again at the expense of assuming, for some 
low Z binaries, sub-primordial helium abundances. Also, one 
may argue that the problem lies in an incorrect Tofi scale 
for low Z stars. Without entering here a detailed discus- 
sion on the robustness of the IRFM scale by Casagrande 
et al. (2006, 2010), suffice here to mention that this is one 
of the hottest scales available (as discussed in the original 
papers and in Sousa et al. 2008). It is hotter by 100 K 
than other IRFM renditions (the now superseded scales by 
Alonso, Arribas & Martinez-Roger 1996 and Ramirez & 
Melendez 2005), and it is comparable within ±50 K to vari- 
ous spectroscopic scales in the metallicity range relevant for 
this work. 

Therefore, any other T^s scale will just worsen the K 
dwarf problem, with real stars even cooler than the models. 
Also, it does not seem plausible that all the (independent) 
scales available are systematically offset by > 200 K, which 
is what is needed to bring the K dwarfs of Casagrande et al. 
(2007) in line with stellar models. 

We estimate the possible extent of the required revi- 
sion of low MS stellar models on the base of homology re- 
lations. First we show that theoretical homology relations 
properly reproduce the response of stellar models to the he- 
lium content: in particular, analyzing the split of the mul- 
tiple MS in cjCen and NGC 2808 by homology relations 
yields the same conclusions as full isochrone analysis (i.e. a 
helium content Y ~ 0.4 for the blue subpopulations) . Then, 
since both isochrones and theoretical homology relations fail 
the interpretation of the nearby low-Z Main Sequence, we 
calibrate empirical homology relations to yield consistently 
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AY/AZ 2 for nearby stars, and inspect the consequences 
for the MS sphts in globular clusters. 

// we entirely impute the failure of the low-MS stel- 
lar models to a wrong response to the helium abundance 
(and correspondingly rc-calibratc the holiurn-scnsitivc term 
of the homology relations), the consequences for globular 
clusters are highly significant: the helium content of the blue 
sub-populations is reduced from 40% to 30%, which is fax 
easier to explain with chemical evolution models (Renzini 
2008; Yi 2009). 

Alternatively, if stellar models for the low MS are as- 
sumed to fail in their metallicity dependence (i.e. we re- 
calibrate only the homology term expressing the metallicity 
dependence of opacity) the consequences for globular clus- 
ters are negligible — as the metallicity differences between 
the subpopulations are minimal or vanishing. 

As the solution for K dwarf models can be intermediate 
between these two extreme assumptions, we suggest that the 
helium rich populations in globular clusters are likely to have 
a helium content Y in between 0.3-0.4; but altogether, there 
is room to decrease their estimated helium content from the 
extreme Y = 0.4 that is the commonly quoted value. 

In summary, the purpose of this exorcise is to draw the 
attention of model-makers to the problem of the HR dia- 
gram of nearby low-Z stars: improvements in this respect 
are potentially important also for the riddle of the helium 
self-enrichment in globular clusters. 

We remind indeed that, while a helium content as high 
as Y ~ 0.35 — 0.4 seems to nicely account for the multi- 
ple Main Sequences and the morphology of the Horizon- 
tal Branch (HB), other observations provide some counter- 
evidence for such a helium rich sub-population. In wCen, 
the location of the RGB bumps of the metal-poor and 
metal-intermediate populations is consistent with a maxi- 
mum Ay < 0.1, less than what is derived from MS and 
HB analysis (SoUima ct al. 2005). SoUima et al. (2006) have 
also identified RR Lyrae stars with metallicity correspond- 
ing to that of the blue MS, but normal helium content — 
which would require the metal-intermediate population to 
be further split into a helium rich and a helium normal sub- 
components. 

NGC 6752 is another cluster suspected to host a he- 
lium enriched population, due to the morphology of its 
HB and the presence of a broadened, possible multiple MS 
(Milone et al. 2010); however, Villanova et al. (2009) did 
not detect large helium abundances in the spectra of blue 
HB stars in the (narrow but crucial) temperature range 
8500 < Tefi < 11500, hot enough to produce helium lines 
but still unaffected by helium sedimentation. Other stud- 
ies of clusters suspected to host helium rich sub-populations 
have not confirmed their presence (e.g. Lee et al. 2009 for 
NGC 1851; Catelan et al. 2009b for M3). While awaiting 
for independent proof of high helium abundances from direct 
spectral measurements or other features in the HR diagram 
(see the recent review by Catelan et al. 2009a), we remark 
that the helium abundances obtained by our proposed revi- 
sion, Y < 0.3, axe consistent with all of the above-mentioned 
constraints. 

The need for sub-primordial helium abundances to fit a 
handful of nearby subdwarfs on the HR diagram was first no- 
ticed by Fernandes et al. (1998) and Lebreton et al. (1999) 
who advocated the inclusion of additional physical processes 



not implemented in standard stellar evolutionary calcula- 
tions. In the following, we discuss a number of possible solu- 
tions to the problem. Broadly speaking, we can classify them 
Eis metallicity-dependent solutions, which shall not concern 
the helium rich populations of globular clusters, and helium- 
dependent solutions, which we expect also to impact the in- 
terpretation of globular clusters. We may add that, consid- 
ering the nice agreement between the multiple MS and the 
HB morpology in uj Cen and NGC 2808, the optimal helium- 
dependent solution should preferably preserve this relation; 
this can be achieved if the new stellar models will not just 
shift the low metallicity ZAMS to lower effective tempera- 
tures, but also affect the luminosity of individual stars and 
accelerate their evolution, which helps to populate the blue 
side of the HB. Casagrande et al. (2007) discussed a number 
of possible solutions (all in the metallicity-dependent cleiss). 
We briefly recall them and discuss a few more here below. 

Mixing length A mixing length parameter decreasing 
with metallicity would make low Z models redder; for the 
isochrones we used, the variation should be from the refer- 
ence solar value a=1.68 to a=1.0 at low Z. While any depen- 
dence of a on metallicity, mass or other physical parameters 
is still highly disputed (Casagrande et al. 2007 and refer- 
ences therein) this solution would have no consequence for 
the helium content of globular clusters, as it is metallicity 
dependent and does not significantly affect the luminosity 
and lifetime of the stellar models. 

Diffusion -|- non— LTE effects Element diffusion affects 

the location of low MS stars in the HR diagram: it makes 
stellar models redder, and moreover it lowers the measured 
surface metallicity with respect to the real intrinsic one. In 
this scenario, low Z subdwarfs look redder than we expect, 
mostly because they are actually more metal rich than we 
measure. Casagrande et al. (2007) discussed this possibil- 
ity resorting to literature models with fully efficient ( "maxi- 
mal" ) difi'usion, and found that it does not completely solve 
the problem, unless it is combined with additional errors on 
the observed metallicity due to non-LTE effects. (This so- 
lution also would have no effect on globular clusters, since 
it is metallicity and age dependent, and the gap in both is 
small for the multiple populations of globular clusters.) 

However, as discussed in that paper, there is evidence 
that both of these mechanisms are not fully efficient in real 
stars, as one would need to solve the K dwarf problem. 
For diffusion in particular, observations of field and cluster 
metal-poor dwarfs point toward inhibited efficiency and ad- 
ditional — yet ad hoc — processes are nowadays invoked 
to contrast diffusion (e.g. Chaboyer et al. 2001; Richard, 
Michaud & Richer 2002, 2005; Korn et al. 2007). 

Boundary conditions also play a role in stellar evolution- 
ary tracks. A solar scaled T — r relation provides a good 
description of the atmospheric structure also in metal-poor 
stars, and may contribute to improve low Z models with re- 
spect to gray atmosphere boundary conditions. However, the 
effect appears to be quite small for MS stars, when each set 
of isochrones consistently adopts a mixing length parameter 
calibrated on the solar model (Vandenberg et al. 2008). 

Besides, boundary conditions are expected to become 
more relevant at lower stellar masses and luminosities (due 
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to the deeper surface convection), while in Casagrande 
et al. (2007) we verified that there is no systematic trend 
of the derived AY/AZ with Mtoi- 

Also this solution relies on a metallicity-dependent ef- 
fect and would not bring substantial revisions to the helium 
content in the multiple MS of GCs. 

Convection One can always consider changing the bound- 
ary of convective regions as a viable working hypothesis. 
(In fact, the very presence of a convective envelope renders 
real K dwarfs not rigorously representable by homology rela- 
tions.) Considering that the convective envelope gets thinner 
at lower metallicities, any change in the convection scheme 
(e.g. extra-mixing, or undershooting) is expected to affect 
solar metallicity objects more, creating a differential effect 
in metallicity that could change the relative location of the 
low MS as a function of Z. Also this solution falls in the 
metallicity-dependent category and is not expected to im- 
pact on the GC analysis. 

Opacity As K dwarfs are largely radiative structures, in 
particular in the regions where nuclear energy is produced, 
opacity is a leading ingredient in determining their luminos- 
ity. An incorrect metallicity dependence of opacity would 
seriously affect the broadening of the low MS. The idea is 
tempting, as the "break" in the estimated AY/AZ occurs 
close to Zo = 0.01, which is the classic divide between free- 
free and bound-free dominated opacity in stellar interiors 
(Cox & Giuli 2004). The parameter P3 in the original ho- 
mology relations is related to this characteristic metallicity 
as P3 = 1/^0 (Section m and Appendix) and we have seen 
that modern isochrones are better described by P3 ~ 670 
rather than 100, effectively reducing Zo to 0.0015 (Section^ 
Pagel & Portinari 1998). However, optimizing empirical ho- 
mology relations on the opacity term, restores a value for 
P3 or Zff/f,f closer to the theoretical one (P3 = 150 or 
Zo ~ 0.007). This suggests that the culprit might be the 
metallicity dependence of opacity (which would then be ir- 
relevant for globular clusters). It is easier to figure errors in 
the bound-free contribution to opacity, rather than in the 
free-free component, at least for the bulk of the star wher 
H and He are completely ionized. Lowering the bound-free 
contribution to opacity would imply a recalibration of solar 
metallicity isochrones, after which high Z (bound-free dom- 
inated) and low Z (free-free dominated) model MS may fall 
closer to each other, thus reducing the broadening. How- 
ever, we wonder if there is much room for profound changes 
in opacity, given the excellent agreement between the ma- 
jor current, independent databases (Opacity Project, Seaton 
2005 and references therein; and OPAL, Iglesias & Rogers 
1996 and references therein) and keeping in mind that im- 
provement in atomic data, input physics etc. most often 
leads to an increase in the opacity, as more and more opacity 
sources are taken into account. Helioseismology has shown 
extant OPAL and OP opacities to be fully adequate for So- 
lar models with the "classic" solar composition; while an in- 
crease in opacity by 10-20% in a suitable temperature range 
has been invoked for the Solar model with new, lower metal- 
licity — possibly in connection with increased neon abun- 
dance (Basu & Antia 2008; Asplund et al. 2009; Serenelli 
et al. 2009; and references therein). How these or other fun- 
damental changes of the Solar model, to recover agreement 



with helioseismological constraints, will impact low Z subd- 
warfs and the relative location of MS of subsolar metallicity 
on the HR diagram, remains to be explored. 

Finally, we notice that the "helium problem" may not be 
limited to K dwarfs: systematic temperature offsets from the 
theoretical Main Sequence at low Z have been highlighted for 
FG dwarfs in the Geneva-Copenhagen survey (Nordstrom 
et al. 2004); and a helium content Y=0. 23-0.24 (i.e. slightly 
sub-cosmological) has been recently suggested for a slighly 
metal-poor F dwarf binary ([Fe/H]=-0.25; Clausen et al. 
2010). If the problem indeed extends to FG dwarfs, clearly 
some of the suggested solutions (for instance, those related 
to diffusion) are not viable. 

All of the solutions suggested above rely, more or less indi- 
rectly, on metallicity dependence; this is the standard way 
we think of stellar models. However, our exercise in this 
paper highlights that it is worth thinking of other possible 
systematics in the stellar models, especially connected to 
the helium abundance, for their interesting impact on the 
multiple MS of globular clusters. 

As mentioned above, opacity has a key role in the struc- 
ture of lower MS stars, so we may wonder whether the issue 
can be the helium dependence of opacity. This can be hardly 
modified in regions of complete ionization; more promising 
are regions where He is partially recombined and contributes 
to the boud-free opacity. 

Other helium-dependent effects may be considered, 
such as diffusion mechanisms that would act differently for 
helium as for metals. (However, preliminary tests on only- 
helium diffusion seem to go in the opposite direction as 
needed for globular clusters, with the helium rich MS getting 
proportionally more red than the helium poor; A. Serenelli, 
priv. comm.) 

Possibly, other of the above mentioned solutions (mix- 
ing length, convection etc.) can be recast and explored in 
terms of helium dependence, rather than metallicity depen- 
dence. Any ideas in this direction are highly desired: hope- 
fully, astrophysicists will be as creative in solving the prob- 
lem of the HR diagram of low metallicity K dwarfs, as they 
have been in tackling the riddle of the extreme helium rich 
populations in globular clusters! 
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APPENDIX A: HOMOLOGY RELATIONS — 
THE BACKGROUND 

An excellent introduction of the homology relations used in 
this paper is given by Fernandes et al. (1996), which we fol- 
low. A family of permanently homologous stars (i.e., stars 
with the same relative mass distribution, that are in hydro- 
dynamical and thermal equilibrium) can be obtained in the 
hypothesis that the equation of state is the perfect gas law. 
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and that opacity and energy generation rate follow laws of 
the kind: 



AAfeoi = -2.5 log A ^ _2.51og i'^^'y ' 



hi = HQ p 1 
e = eo p 1 

(Cox & Giuli 1968). Other underlying assumptions are that 
the energy transport occurs via radiative diffusion over most 
of the interior (in particular, over the energy generation re- 
gion, as is the case for the low MS) and that the radial profile 
of chemical composition differs, in different stars, just by a 
scale factor. The latter assumption holds in particular for 
stars on the ZAMS, with internally uniform chemical com- 
position. 

For low-mass MS stars, opacity is dominated by bound- 
free and free-free processes and well approximated by 
Kramer's law: k — K,apT~^'^\ while the energy generation 
rate via pp chain approximately follows epp = eopT"^. With 
these dependencies on p and T, homology transformations 
yield the following relation for the luminosity and effective 
temperature of the star 



T -0.077 7.769 -1.077 ,«-5. 462 

L = en p- Kn M 



Ton = e 



g /J, n,Q 

0.096 ,2.211 .,-0.346^1.327 



as a function of mass M and molecular weight p. Eliminating 
mass, one derives: 



r 0.318 -1.332 0.347 rrn4. 116 

— eo M '^o -'off 



(Al) 



The last factor indicates the slope of the ZAMS for a given 
chemical composition: about 10 mag per dex in log(reff ) — 
flatter than, but not so far from, the slope obtained with de- 
tailed stellar structure computations (17 mag per dex, Sec- 
tion [2]). Here we are mostly interested in the broadening of 
the ZAMS in luminosity as a function of chemical composi- 
tion, at a fixed Tcff ; such broadening is expressed by the first 
three factors. Since H-burning occurs via the p-p chain: 

eo oc A 

The molecular weight p, throughout most of the star where 
complete ionization applies, is given by: 

p-^ ^2X + ^Y+lz ^ i(3 + 5X-Z) 
4 2 4 

For bound-free and free-free opacity, approximately: 
no,bf « 7.40 X 10^'' X -^Zil + X), Hiojf ~ 3.76X 10^^(1 + X) 
(Cox & Giuh 2004), so that 

Ko = fio,bf + no.ff oc (1 + X){100Z + 1) 

Notice that Zq — is where opacity drifts from free-free 
dominated {Z < 0.01) to bound- free dominated (Z > 0.01). 
Altogether, we can write the homology relation (|A1|) as: 

L = F{eo,p, Ko) ToB^^ 

with 

T a X''-^^^ (3 + 5X - Zf-''^^ (1 -f Xf'^^lOQZ + 1)°-^*^ 

The first factor in the T function is related to the energy 
generation rate, the second factor to the molecular weight, 
and the last two factors to opacity. It is then straightforward 
to derive Eq.[l]for the broadening of the ZAMS in luminos- 
ity, with respect to a reference composition (Xr,Yr, Zr): 



Al Polynomial approximations 

Pagel & Portinari (1998) suggested an approximation of the 
homology relations based on Taylor series expansion of the 
terms in log(l — x): 



log(l - a;) 



1 

In 10 n 

n = l 



Ixl < 1 



Expansion to first order of Eq. [5] yields a linearized formula 
similar to eq. 5 of Pagel & Portinari: 



AlogTeff(i) = _o.051 log ^ 



lOOZr + l 



I I 0.0406 I 0.02215 _, 0.4256 



' l+Jfr ' 3+5X,,-Z; 

0.08512 



^5{Z~Zr) (A2) 



:{Z^Zr) 

(where the last term is negligible). However this linearized 
approximation holds only for \5{Z — Zr)\ « 1 and does not 
recover very well the original homology relations for > 5 
(see Fig. Al). 

One can then resort to further terms in the expansion. 
The second order terms yield: 

A log Toff = A log Toff 

, [ 0.0203 I 0.01 



107 I 1.064 X'-^t'y V 

XTP' + (3+5X,.-Z,.)2 I " - 



I 0.4256 a. I '7 7^ 

+ (3+5X,-Z,)2 <^\^ ~ ^r) 



+ 



0.04256 



(3+5Xr-Z^)- 



{Z - Zrf 



where the last two terms are negligible. The third order ex- 
pansion terms yield: 

AlogTeff'^' = Al0greff'2^ + 



01354 I 0.00738 , 3.547 



'W 



(1 + XrF (3+5Xr-Zr)^ 



] 5\Z~Zr) 



0.2128 x2 I- 7 7 \3 
(3+5X,-Z,)a \^ - 

(3+5X^-2^)7 — Zr)^ 



(3+5X,.-Z,.)'3' 
0.0284 



where the last three terms are negligible. In general, neglect- 
ing the smaller terms of the kind 5'"(Z — Zr)^ with m < n, 
we can write the following polynomial approximation to ho- 
mology relations: 



AlogTeff- -0.051 log (^f|±^) 



. Y^oo [ 0.0406 I 0.02215 _i_ 
~rZ^„=l [ x;> + (1+X^)" 



(3+5X 



0.08512 1 i5"(Z-Zr) 
5Xr-Zr)"J n 



(A3) 



The corresponding parametric fitting formula for the 
isochrones would be of the kind: 

AlogTeff --ailogf "'f + M +f2b„5"{Z-Zrr (A4) 

V a2Zr + 1 / — ■ 

^ ^ n — 1 

For n-th order series expansion, this fitting formula has n + 
2 free parameters. Pagel & Portinari (1998) indeed used a 
linearized approximation {n — 1) similar to Eq. IA2I as a 
guideline, and had 3 parameters in their fitting formula for 
isochrones. 



14 Portinari, Casagrande and Flynn 




0.005 0.01 0.015 0.02 0.01 0.02 0.03 

Z Z 



Figure Al. Homology relations for ranging from to 20 (with a pace of 2) as indicated. Coloured dashed lines represent the 
approximate relations of Eg. IA3l for increasing order of the terms of the expansion. The true relations are well reproduced up to = 20 

for n > 2 — at least within a metallicity range AZ ~ 0.02. The dotted line shows our simplified form of the homology relations (Eq. [Sjl . 



Fig. Al shows that expansion to order n > 2 is needed 
to yield a good fit up to AY/ AZ ~ 10 (at least in the range 
Q ^ Z < Zq) needed for the study of nearby K dwarfs of 
low Z. Our simplified form of the original homology rela- 
tions (Eq. [S]) performs as well as the n = 3 order of series 
expansion, with just 3 free parameters in the corresponding 
fitting formula. 

For extremely high AK/AZ, as applies to the case of 
cjCen, it is worth using the original formulation of the ho- 
mology relations (Eq. (2]), which indeed we showed in Sec- 
tion [3] to yield excellent results when applied to the multiple 
MS of wCen amd NGC 2808. Notice, in fact, that for very 
high Ay/ AZ the basic condition for Taylor series expansion, 
\x\ = \5{Z — Zr)\ < 1, may no longer be fulfilled. 



APPENDIX B: UPDATED HIPPARCOS 
PARALLAXES 

For simplicity this paper relies on the helium abundances 
determined by Casagrande et al. (2007); in the meantime, 
the HR diagram of nearby stars has been revised with up- 
dated Hipparcos parallaxes (van Leewven 2007). In Fig. Bl 
we show that the effect of the revised parallaxes is minimal: 
the sample only slightly changes (3 stars are now rejected 
by the parallax limit 6%, while 5 new stars are included, 



for a total of 88 objects) and the trend toward significantly 
sub-cosmological Y values is confirmed. Nor it appears to 
be due to low metallicity stars being systematically more 
distant than the solar metallicity ones (with systematically 
intrinsic bluer colours due to neglected reddening) : though a 
slight distance-metallicity trend is present as expected, most 
of the stars lie within 30 pc for all metallicities (Fig. B2). 
Also, in Fig. [1] we show that the reddening/extinction cor- 
rection needed to reconcile metal-poor stars with isochrones 
would correspond to E{B — V)=Q.1Q, which is certainly too 
extreme for so nearby stars. 

Other possible systematics have been carefully checked 
and excluded in Casagrande et al. (2007). We do not deem 
necessary here to repeat the detailed Monte-Carlo simu- 
lations performed in 2007 to estimate realistic errorbars 
(shown in Fig. [5} , since it suffices here to demonstrate that 
the new parallaxes by no means solve the riddle of sub- 
cosmological helium abundances. 

Finally, we have also verified that the updated temper- 
ature scale by Casagrande et al. (2010) also does not signif- 
icantly change the Y(Z) plot. 
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Figure Bl. Helium plot for the sample by Casagrande et al. 
(2007) (stars) versus the update sample with the revised paral- 
laxes of van Leewven (2007) (squares) 
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Figure B2. Distances and metallicities of sample stars from 
Casagrande et al. (2007) (top panet) and for the updated sample 
with the revised parallaxes of van Leewven (2007) (bottom panel) 



